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Abstract. We consider the impact of thermal inflation - a short, secondary period 
of inflation that can arise in supersymmetric scenarios - on the stochastic gravitational 
wave background. We show that while the primordial inflationary gravitational wave 
background is essentially unchanged at CMB scales, it is massively diluted at solar 
system scales and would be unobservable by a BBO style experiment. Conversely, 
bubble collisions at the end of thermal inflation can generate a new stochastic 
background. We calculate the likely properties of the bubbles created during this 
phase transition, and show that the expected amplitude and frequency of this signal 
would fall within the BBO range. 
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1. Introduction 

In the coming decades, it is very likely that gravitational waves will be observed by 
direct-detection experiments. These include the presently operating LIGO and VIRGO 
observatories, and proposed experiments such as LISA [1], BBO [2] and DECIGO [3]. 
It is likely that the first sources to be detected will be locahzed, transient events, but 
many cosmological processes yield a stochastic background of gravitational waves. These 
backgrounds are a key target for the nascent field of gravitational astronomy, since they 
probe the primordial universe in ways that are presently impossible. 

The most widely discussed cosmological signal is the inflationary gravitational wave 
background, generated by quantum fluctuations of spacetime. The amplitude of this 
signal decreases with the inflationary energy scale, rendering it effectively invisible in 
models where inflation occurs well below the GUT scale. A stochastic background 
of gravitational waves can also be generated whenever the universe undergoes a 
phase transition. Unlike the inflationary background, which is scale invariant, these 
backgrounds will be peaked at a wavelength correlated with the characteristic scale of 
the transition. In these scenarios, the gravitational waves are sourced by the non-zero 
quadrupole (and higher) moments of the matter distribution [4]. 

Recently, several groups have investigated the possibility that gravitational waves 
are produced during the preheating phase at the end of inflation [5, 6, 7, 8, 9, 10, 11]. The 
momentum modes of fields coupled to the inflaton undergo resonant pumping, which 
renders the universe highly inhomogeneous, and sources the emission of gravitational 
radiation. Intriguingly, this signal is easier to detect in models where the inflationary 
scale is substantially below the GUT scale, as first argued in [6], and confirmed in [7]. 
Moreover, it probes the mechanism that ends inflation, which is normally assumed to 
be inaccessible to direct observational tests. 

In addition, a second (and not entirely unrelated) mechanism could generate a 
significant background - namely the collision of percolating bubbles formed during a 
first order phase transition [12, 13, 14, 15, 16, 17]. As we will see, it has the same 
scaling properties as the preheating signal - namely the amplitude need not depend 
on the energy scale, but its (present-day) wavelength decreases linearly with this scale. 
Assuming an energy scale somewhere between the electroweak and GUT scales, the 
Hubble scale will he between 10~^ eV and 10^^ GeV. We will see that if the characteristic 
bubble radius is within a few orders of magnitude of the Hubble scale at which they 
form, the resulting gravitational wave background will have a present-day frequency in 
the range 10"''^ to 10^ Hz [6, 16]. 

To predict the primordial background, one must know the transfer function, which 
is determined by the effective equation of state of the universe between the epoch when 
the signal is generated and the present day. A further comphcation is that there is 
no guarantee that inflation happens only once in the early universe. One may imagine 
primary and secondary periods of inflation - where the primary period is responsible 
for solving the usual cosmological initial conditions problems and laying down the 
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perturbations that seed large scale structure, and the secondary period lasts for < 0(10) 
e-foldings, and occurs at some energy scale below that of the first inflationary period 
and above (or at) the scale at which baryogenesis takes place. If signiflcantly more than 
(9(10) e-folds occurs in the secondary phase, it will erase all relics of the previous phase, 
returning us to what is effectively a single phase model. In addition to modifying the 
mapping between scales in the present epoch and the moment these (comoving) scales 
leave the horizon during the primary phase, a secondary phase of inflation also erases 
the primordial inflationary background of gravitational waves at wavelengths which have 
reentered the horizon before the secondary phase begins. This leads to a situation where 
a high scale primary phase would generate a gravitational wave background visible 
in the B-mode of the CMB polarization, but completely invisible to direct detection 
experiments operating at solar system scales, such as BBO. 

A cosmological model with multiple inflationary phases may appear baroque. 
However thermal inflation [18, 19] is a secondary period of inflation that arises in 
super symmetric theories when some (almost) flat direction with negative mass-squared 
at the origin - which we will call the flaton - acquires a thermal contribution to its 
effective potential, giving it a positive effective mass. The flaton is thus pinned at the 
origin, and the flnite potential energy at this point drives a secondary period of inflation. 
Once thermal inflation sets in, the temperature falls rapidly and when it drops below 
the flaton's mass scale (usually assumed to be a typical soft supersymmctry breaking 
mass scale ~ TcV) the flaton rolls aways from the origin, and thermal inflation naturally 
shuts off. The total amount of inflation is typically small, about 10 e-folds or less. 

Thermal inflation requires only a flat direction with negative mass-squared at 
the origin, both of which are natural in supersymmetric models. Prom the model 
building perspective, its importance lies in its ability to dilute the moduli and gravitinos, 
which are generically produced after the primary period of inflation in supersymmetric 
theories. These particles either decay, wrecking the successful predictions of Big Bang 
Nucleosynthesis, or provide a dark matter density that massively overdoses the present 
universe [20, 21, 22, 23, 24]. These particles are produced before thermal inflation 
begins and are too weakly coupled to be produced after it ends, so their contribution is 
thus diluted to safe levels. Thermal inflation is perhaps the best-motivated solution to 
the moduli problem, which is endemic to supersymmetric models of the early universe, 
including many of those derived from string theory. Thermal inflation has a relatively 
low reheating scale (typically O{10) GeV) but it can give rise to baryogenesis [25, 26]. 

Unfortunately, as we will show in Section 4, thermal inflation wipes out any 
primordial inflationary gravitational wave signal in the range of frequencies accessible 
to BBO |. Thus, even if the primary phase of inflation occurs at a high enough scale 
to produce a detectable tensor background, thermal inflation erases it at short scales 
[27]. If this was the end of the story, one might regard it as an unfortunate side-effect 

I Indeed it also wipes out the primordial density spectrum at the same frequencies, although the 
scales affected are not important cosmologically and hence not constrained by large scale structure 
observations. 
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of solving the moduli problem, but there is a final piece to the puzzle: thermal inflation 
ends via bubble nucleation, and the collisions between these bubbles as they percolate 
produces a significant stochastic gravitational wave background. We will see that this 
process may generate i^cw^^l/ ~ 1 Hz) ~ 10^^'', and is potentially detectable by BBO. 

Thermal infiation is a complicated scenario, at least when compared to simple 
toy models of inflation. However, it relies on generic ingredients found in many 
supersymmetric or stringy models of the early universe; it is the intrinsic richness of 
these models that generates their complicated phenomenology. In particular, the moduli 
problem is robust enough to have survived 25 years. Conversely, thermal inflation does 
not require any exotic physics, unnatural parameter values or initial conditions, and 
the complexity of this scenario may in fact be typical of realistic models of the early 
universe. Finally, while the primordial gravitational background is erased, it would need 
to have been near the top of its plausible range to have been observable even without 
dilution, and furthermore in supersymmetric and stringy models the bottom of this 
range is arguably more natural [28, 29, 30]. The bubble collision signature is generic to 
thermal inflation and effectively expands the hst of the models probed by BBO. 

In Section 2 we summarize thermal inflation, and describe the process of bubble 
nucleation that occurs as it ends. In Section 3 we describe the production of gravitational 
waves by bubble collisions and estimate the expected spectrum. In Section 4 we compute 
the suppression to the primordial inflationary background by a secondary period of 
inflation, combine this with the likely spectrum from bubble collisions, and assess the 
prospects for detecting this signal in future observatories, the likelihood it can be 
distinguished from other stochastic astrophysical backgrounds, and the likely tensor 
signal from the primary inflationary phase, which would be undiluted at CMB scales, 
or via pulsar timing experiments [31]. Finally we sum up in Section 5. 

2. Thermal inflation 

2.1. Motivation 

The main motivation for thermal inflation [18, 19, 32, 33, 34, 35, 36, 37, 38] is to solve 
the moduli problem [20, 21, 22], though it also solves the gravitino problem [23, 24] and 
provides a mechanism for baryogenesis [25, 26, 39, 40, 41, 42, 43, 44]. Irrespective of 
these uses, thermal inflation is sufficiently natural that it might occur anyway. 

Moduli are scalar fields with Planckian vacuum expectation values, and hence 
gravitational strength interactions. Their potential arises due to supersymmetry 
breaking, and assuming supersymmetry breaking is transmitted to the observable sector 
via gravitational strength interactions, moduli have vacuum masses of order the soft 
supersymmetry breaking scale in the observable sector 

m„,oduii ~ ~ 10" to 10^ GeV . (1) 

However, in the early universe, the finite energy density breaks supersymmetry. When 
H ^ rUg this supersymmetry breaking dominates over the vacuum supersymmetry 
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breaking and hence determines the moduh potentiaL When H drops below nis the 
moduh potential reduces to its vacuum form, but with the moduli typically displaced 
by a Planckian distance. The moduli then start oscillating with Planckian amplitude 
and immediately dominate the energy density of the universe. Thanks to the relatively 
low moduli mass and their very weak interactions, these oscillations persist beyond 
nucleosynthesis with disasterous consequences [20, 21, 22]. 

Inflation is typically invoked to rid the universe of unwanted relics, but one would 
expect an inflaton to have a mass > rus and so primordial inflation to occur at scales 
H > rUs- The moduli are generated at if ~ mg, and to a lesser extent by any phase 
transition at H < rus- Thus one wants inflation at if -C mg to dilute the moduli, but 
it is very difficult to realize primordial inflation at these scales. On the other hand, 
thermal inflation [18, 19] automatically occurs at if -C mg. For a thermal inflation scale 

10^ to lO^GeV (2) 

thermal inflation provides enough dilution to rid the universe of moduli, but has a low 
enough scale not to regenerate them afterwards. Furthermore, as thermal inflation lasts 
about 10 e-folds or less, it does not destroy the primordial perturbations needed for 
structure formation. While thermal inflation provides a very natural solution to the 
moduli problem, it is incompatible with most baryogcncsis scenarios since it will dilute 
baryons produced before it begins, and the reheat temperature after thermal inflation 
is very low (0(10 GeV)). Fortunately, a baryon asymmetry is almost automatically 
generated at the end of thermal inflation [25, 26], solving this difficulty. 



2.2. Particle physics model 

Like the Standard Model Higgs field, the flaton has a negative mass-squared at the 
origin 

with m<^ ~ mg, the soft supersymmetry breaking mass scale. Unlike the Standard 
Model Higgs field, but like many scalar field directions in the Minimal Supersymmetric 
Standard Model (MSSM) §, the flaton does not have a stabilizing 0^ term. Instead, non- 
renormalizable higher order terms stabilize its potential at a large field value 0o ^ ^tt.^, 
as illustrated in Figure 1. The potential energy at the origin 

Fti ~ mlcg (4) 

is tuned to give zero vacuum energy density. For thermal inflation to solve the moduli 
problem, we require [19] 

KrY^ ~ 10^ to lO^GeV , (5) 

§ For example, if the MSSM constraints m^^ + m|^^ + 2\n\'^ > 2 \Bii\ and (rnjj^ + — m\){m\ + 
rrvfj + l/ip) > |-B/i|^ were not satisfied, then the MSSM would have a flaton. 
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Figure 1. Thermal inflation occurs wlien a flaton (f> is held at the origin by its flnite 
temperature potential and Vti dominates the energy density. 



corresponding to 

00 ~ 10^° to 10^2 GeV . (6) 

The flaton is a Standard Model gauge singlet, since it would otherwise break the 
Standard Model gauge group at the scale 0o; but has unsuppressed Yukawa couplings to 
non-MSSM fields \E' and ^, which acquire masses Mij, ~ 0o once the fiaton reaches the 
minimum of its potential. \1' and ^ are constrained by gauge coupling unification, and 
at least some of their Yukawa couplings to the fiaton should be strong enough to drive 
the fiaton's mass-squared negative at low energies, in the same way as the top Yukawa 
coupling drives the Higgs mass-squared negative at low energies. Such a strong Yukawa 
coupling will be attracted to its renormalization group infrared fixed point, allowing one 
to relate its value at low energies to the strong gauge coupling and the choice of and 
^ representations. See Appendix C for more details and numerical estimates for this 
Yukawa coupling. 

The fiaton's Yukawa terms couple it strongly to the thermal bath, so its finite 
temperature effective potential is 

V(^,T).V.,J-f:t\^^^-<)^^- f ^<<^ (7) 
1 -^m^(p^ + . . . for > T 

where the values of A and B depend on the details of the thermal bath and the fiaton's 
coupling to the thermal bath, as we will see in Appendix B. The potential in Eq. (7) 
has a global minimum at = 0o for T <^ 0o and a local minimum at the origin for 

T > = ^ (8) 

V B 



where Tc2 is the critical temperature for the would be second order phase transition. 
However, we will find that thermal infiation always ends slightly before the temperature 
drops far enough for the quadratic term to change sign. The field is thus trapped 
at the origin and the phase transition is first order, leading to bubble formation at a 
temperature Td > Tc2. Since Td ~ Tc2, when the difference between the temperatures 
is unimportant we will simply denote them by Tc. 
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Figure 2. History of the universe witli tliermal inflation. The fractional density 
Q.X = Px/p of moduli — , potential — , flaton — and radiation — is plotted against 
the number of e-folds of expansion In a. 



2.3. Cosmological history 

The overall history of a universe which undegoes thermal inflation is summarized in 
Figure 2. The primordial inflationary phase, which lays down the perturbation spectra, 
ends at some time t < t^. In supersymmetric theories, the moduli are produced at the 
time t = tra when H ~ mg. At this time the moduli potential switches from its finite 
energy density form to its vacuum form, leaving the moduli oscillating with Planckian 
amplitude and PmoduU ~ Prad- However, the moduli quickly dominate, as Pmoduii oc 
and Prad oc a~^. 

Thermal inflation begins at tb, at which point Pmoduii ~ Vti and 

Pmoduli _ f _ f TT^s^Pl\ ^g^ 



Prad (^la 



where Mpi = l/VS^ = 2.4 x 10^^ GeV, and 



n-(^Y'vil'r.(jS^Y\ (10) 

V Pmoduli/ Xm^MpiJ 



Thermal inflation ends at tr when 



c 

T = Tcr^mj,. (11) 
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During thermal inflation, the temperature of the thermal bath drops as T oc a~^, 
assuming the effective number of degrees of freedom does not change significantly. 
Therefore, from Eqs. (10) and (11), thermal inflation lasts for less than 10 e-folds ||, 
since 



A^Ti - In 



6 



In 



0, 



6 to 9. 



(12) 



After thermal inflation, the universe is dominated by the flaton and oc a ^ ^. Finally, 
the flaton decays at t — t^ leaving a radiation dominated universe at temperature 

Td ~ 1 to 10^ GeV . (13) 

At this point, the stage is set for the standard hot Big Bang. 



2.4- The first order phase transition at the end of thermal inflation 

The origin is always only a local minimum of the flaton's finite temperature effective 
potential, all the way down to the temperature Tc2, so the phase transition that ends 
thermal infiation will be first order. We therefore apply the thermal tunnelling formalism 
of Refs. [45, 46, 47]. After nucleation, there is a very large pressure difference across 
the bubble wall due to the deep vacuum of the fiaton field, hence we expect the bubble 
wall propagation to proceed as a detonation [48], and the walls rapidly accelerate to 
relativistic speeds. Finally, the bubbles will percolate, ending thermal infiation. 
The bubble nucleation rate is 

SsiT)- 



r{T) - T^exp 



T 



(14) 



where Ss{T) is the energy of a bubble. The temperature Td at which the dominant 
bubbles nucleate is given by 

r(TeO'' 



r(rci) 



r(T, 



cl 



(15) 



see Appendix A. The timescale of the transition, from nucleation to percolation, and 
hence the characteristic size of the bubbles at percolation, is given by the important 
parameter (3^^ 

T 



P 



where 



Hr 



HV 



Tel 



1 d 
Hdi 



nn d 

T — 

dT 



S3 
T 



(16) 



TI 



3M2j 



(17) 



is the Hubble parameter at the time of nucleation and percolation. The exact value 
of (3 depends on the strength of the phase transition, which in turn depends crucially 

II This estimate could be increased by about 5 e-folds in the baryogenesis scenario of Refs. [25, 26]. 
^ Preheating may alter this equation of state slightly. 
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on the details of the thermal bath which holds the flaton at the top of its potential. 
In Appendix D, we use the constraints on the flaton's Yukawa couplings, obtained in 
Appendix C, to derive estimates for (3 /He, finding 

4- ~ 10^ to 10^ . (18) 
He 

3. Gravitational wave production during bubble percolation 

Gravitational waves are generated when the bubble walls collide [14], and there may be 
an additional contribution from turbulence [16]. The ratio 

(3 _ t 



HF 



(19) 

determines both the peak frequency and amplitude of the gravitational wave spectrum 
generated by bubble percolation. Our result for P/Hc, Eq. (18), tells us that the bubble 
radii at percolation are much smaller than the Hubble radius. This justifies the use of 

the gravitational wave spectrum derived numerically in Ref. [16] which assumed that 
the gravitational waves arc generated in an effectively ffat background. Physically, 
(3 /He determines the peak frequency because the size of the bubbles determines the 
characteristic lengthscale of the quadrupole moment that generates the gravitational 
waves, and hence the frequency at which we observe them today. 

The peak amplitude is a httle more subtle. The total energy radiated in gravity 
waves, per bubble, per frequency interval dcu and solid angle fl is [4] 
dE 

= 2Gu'Kj,im{k)T*^{k,u;)Ti,^{k,u;) (20) 

dual I 

where Aj^ is the projection tensor that picks up the traceless and transverse component 
of the source stress tensor, with k = k/|/c| 

1.... 1 1.. 1 

■^ij,lm — ^ij ^jm 2 kj k^ + — k j k j k; km — 5ij 5im'\"^ ^ij k; IS-m -\- — Simk.ik.j .(2,V) 

Tjj(k, u) is the fouricr transform of the source stress tensor for each nucleating spherical 
bubble with size i? [16] 

Tij{k,u;) ^ — f dt e^^ f dVt R^p^n±i±j. (22) 
Stt Jq J 

K < 1 is the efficiency factor for the conversion of the vacuum energy into the kinetic 
energy of the bubbles. This imply that, per bubble, dEgy^/dw oc {R^KpyY- Notice 
that it is a product of the vacuum energy density which can be traced back to the 
fact that the quadrupole moment of a gravitational source is proportional to the source 
energy density. Meanwhile, the total energy is simply E,^ oc R^Pv, and substituting the 
lengthscale of the relativistic bubble R (x and integrating over the peak production 
frequency lu oc P we get 

^nawocK'f^y. (23) 



-E'total \ P 
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recalling that H'^ oc p^. This relation was first derived semi-analytically and checked 
numerically in references [15, 14]. In particular, the total energy in gravitational waves 
generated by thin-wall bubble collisions is well-fitted by the heuristic formula [16] 

where a is the ratio of vacuum energy to thermal energy before the phase transition 
occurs 

« = ^^-10-, (25) 

V is the wall velocity, and as discussed in Section 2.4 we expect w ~ 1, and since the 
thermal energy is negligible we expect k ~ 1. 

In order to assess the observabihty of this signal by a given detector we need the 
spectrum 

OGw(t^) = -T\ ■ (26) 

This spectrum has been numerically computed for transitions of the sort we are 
concerned with here [15, 16]. In particular, looking at Figure 7 of Ref. [16] with v — 1, 
the peak value for JIqw is 

Qowioo /3) = 0.03 X — (27) 

In other words, the peak power is at a characteristic frequency u = j3. This is then 
rescaled by the subsequent cosmological expansion, so the frequency today is 



/«o) = i^ if - (28) 




0-7 a. r-^ -^ Rfe^ . (29) 



The last factor is unity if the universe is dominated by non-relativistic fiaton matter after 
thermal inflation, but bubble collisions and other preheating effects would be expected 
to heat up the flaton matter at least initially, making this factor greater than unity. 
However, these effects are very difficult to estimate so we assume the most conservative 
value of unity. Assuming a more favorable set of thermal inflation parameters for the 
remaining factors {f3/H^ = 10^, V^^^ = lO^GeV, = lO^GeV), the peak frequency 
today is about 0.7 Hz, well within the spectral range of future space-based gravitational 
wave detectors such as LISA [1] and BBO [2] . Similarly, we can map the power spectrum 
into its present day form 



fiGw(/p.ak,«0)ft' = f!Gw(/c,ic)AM- TT (3") 




8xio-'M^I liT^r] (tt^I 1^1 m 

Pd «d 
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With the same parameters as above, we find a peak gravitational wave power of 
^Gw(/) h'^ 5 X 10~^®, which is close to BBO's design sensitivity, and well below 
that of LISA. 

This amplitude is lower than the i^Gw(/) h"^ ~ 10^^^ seen after preheating or 
some turbulent processes [5, 7, 6, 8, 9, 10, 54, 55]. The bubbles are initially small, 
since (3/Hc ^ 1000, which leads to a shght suppression of power, via Eq. (31). More 
importantly the period of flaton matter domination following bubble nucleation dilutes 
the power by a redshift factor of Oc/od ~ 10~^, relative to what would be seen if the 
universe became radiation dominated immediately after the bubble collisions. On the 
other hand, we have not considered gravitational waves sourced by turbulent mixing well 
after the bubble collisions, which is a potential source of additional power - especially as 
these gravitational waves would be less diluted by the fiaton dominated epoch. Eq. (31) 
is thus a lower hound for the gravitational wave power. Further, much of the energy 
of the universe is contained in the walls, which are moving at relativistic speeds. This 
should help to make the initial preheating of the flaton very efficient, and the details 
of this will need to be carefully studied in order to determine the gravitational wave 
spectrum generated in the post-collision universe. 

In addition, Caprini, Durrer and Servant [17] have recently given an elegant analytic 
discussion of gravitational wave production during bubble collisions. The amplitude and 
shape they obtain for the resulting spectrum is broadly similar to the numerical results 
of [16] (note that [17] assumes there is no matter-dominated phase after percolation), 
adding support to our estimate. They also note that bubble collisions can in fact be 
a sub-dominant source of gravitational waves following a phase transition, and a full 
analysis of this signal will require further study. 

4. Forecasts for Observations 

If a CMB B-mode is sourced by a primordial spectrum of gravitational waves, the 
properties of this signal do not depend strongly on the details post-infiationary physics, 
since the relevant modes arc outside the horizon until after matter-radiation equality. 
Conversely, direct detection experiments are sensitive to modes with much smaller 
wavelengths, and their current amplitude can depend strongly on the expansion history 
of the post- inflationary universe [61, 62]. In particular, since gravitational waves scale 
like radiation with pow oc 1/a^, any phase in which the universe is dominated by 
a component whose density scales away less rapidly than radiation will suppress the 
gravitational background on sub-horizon scales, as first discussed by [27]. Looking at 
the cosmological history implied by thermal infiation, we find three distinct epochs 
during which the short wavelength primordial gravitational wave background can be 
diluted: moduli matter domination, thermal inflation itself, and flaton domination after 
thermal inflation. 

In Figure 3 we show the history of flve representative modes. All modes reenter the 
horizon after the primary phase of inflation, with an amplitude determined by standard 
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Figure 3. Wc illustrate the evolution of five representative gravitational wave modes 
(dotted lines) and the Hubble radius (solid line) in a universe with thermal inflation. 
Physical wavelength A is plotted against the scale factor a. The axes are natural 
logarithm scales normalized so that the scale factor and the Hubble radius are unity at 
the end of thermal inflation. For (3/Hc ~ 10^, the bubbles are about times smaller 
than the horizon. The violet patch indicates the times and scales at which bubble 
collisions generate gravitational waves. The wavelength of Mode 1 is large enough to 
keep it outside the horizon until well after thermal inflation, and it is thus unaffected 
by the turmoil within the horizon. Mode 2 enters the horizon during moduli matter 
domination, and its amplitude is diluted during this phase and during thermal inflation, 
but it freezes out when it re-exits the horizon. Mode 3 resembles Mode 2 except that 
it reenters the horizon during flaton matter domination, and is thus suppressed a 
second time. Mode 4 remains within the horizon during thermal inflation, and is this 
maximally suppressed. Conversely Mode 5 also remains within the horizon, but is 
short enough to be sourced by the bubble collisions. 



inflationary perturbation theory [56], QQ^{k) oc H^^^. During moduli matter domination 
a mode inside the horizon will be diluted by a factor of a~^. As soon as thermal inflation 
kicks in the suppression factor is a~^, dropping back to zero if the mode re-exits the 
horizon and re-freezes. Modes inside the horizon are further diluted by a factor of 
during the fiaton dominated phase, but this extra suppression is largely academic, other 
than for modes which are sourced by the bubble collisions. A schematic plot of the final 
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Figure 4. The expected gravitational spectrum for a universe witli a thermal 
inflationary phase, with projected sensitivities of future planned experiments. The 
black line shows the upper bound on the primordial gravitational wave spectrum from 
CMB data at the time of writing. For reasonable parameter values, thermal inflation 
leads to a suppression of the primordial signal for / > 10^^ Hz, rendering it invisible to 
direct detection experiments. The cyan line indicates the spectrum generated by bubble 
collisions at the end of thermal inflation, using the fiducial values of P/H « 1000, 
Vr^('^ ~ 10^ GeV and Td ~ 10^ GeV. The magenta hue is the expected GW spectrum 
from parameteric resonance [59, 60] after modular inflation [28, 29, 30]. It gets wiped 
out by thermal inflation. 



spectrum is given in Figure 4. Note that the primordial gravitational wave spectrum is 
unsuppressed on long scales, so the constraint on the inflationary scale implied by the 
upper bound on the CMB B-mode is still valid [57]. 

For a fixed value of ^^gw(/)) gravitational waves are easier to detect at long 
wavelengths, as the required strain-sensitivity scales as [58]. Via Eq. (29), we see that 

1/3 

the peak frequency scales as He , and a lower thermal infiation scale leads to a redder 
peak, and thus a more detectable signal. Further, for fixed Td, lowering the thermal 
infiation scale shortens the fiaton dominated period, which, from Eq. (31), increases 
the present day amplitude by a factor Hc^^^. On the other hand, if is increased, 
the amplitude is boosted by a factor of T^^^, while the peak frequency grows as T^^. 
Thus, low scale thermal infiation with a short period of fiaton matter domination is 
the optimal prescription for generating an observable spectrum of gravitational waves 
during the bubble nucleation phase. In terms of the thermal infiation parameter space, 
this situation is equivalent to having a small value of 0o- 

For generic parameter values, this signal will peak at around frequencies of 1 Hz, as 
shown in Figure 4. This corresponds to a wavelength on the order of 10^ km, which is 
characteristic of BBO or DECIGO style space-based proposals [63, 61]. Moreover, the 
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strong spectral dependence of the bubble collision signal will differentiate it from that 
the featureless primordial spectrum. 

5. Discussion 

This analysis draws on a number of apparently disparate topics in theoretical cosmology. 
Firstly, it looks carefully at thermal inflation, considering both its expansion history 
and the properties of the bubbles created by the phase transition that marks its end. 
Historically, the primary motivation for considering thermal inflation is that it provides a 
natural solution to the moduli problem, which is endemic to supersymmetric models and 
arises from the production of a condensate of long-lived moduli particles which either 
disrupt nucleosynthesis as they decay or massively overdose the universe. Thermal 
inflation dilutes these unwanted particles to acceptable levels. However, we show that 
thermal inflation also dilutes the primordial gravitational wave background at small 
scales, rendering it effectively unobservable by direct detection experiments such as 
BBO, but preserving any signal that is present in the CMB B-mode. The detailed 
dependence of the BBO signal on the effective equation of state in the primordial 
universe has been discussed elsewhere (e.g. [61]) but the dilution factor implied by 
thermal inflation would effectively erase any primordial signal. 

Conversely, thermal inflation ends via a flrst order phase transition, so the 
percolation of the resulting bubbles can generate a new background of gravitational 
waves. This spectrum is conflned to a few decades in frequency - unlike the scale- 
free inflationary background - and is potentially detectable by BBO style experiments. 
While the detailed mechanism is different, this signal has much in common with the 
gravitational wave background that can be generated during preheating or parametric 
resonance [5, 6, 7, 8, 9, 10, 11]. Further study will be needed to see if the characteristic 
spectra can be easily distinguished from one another. 

The analysis here only considers the background generated by the bubble collisions 
themselves - a period of turbulent mixing after the bubble collisions would also source 
gravitational waves, and this signal has not yet been computed. It has the potential to 
signiflcantly enhance the gravitational wave signal, since the bubble nucleation phase 
is followed by a period of flaton domination, mimicking a matter dominated universe, 
and diluting the background produced by colliding bubbles. If the turbulence lasts 
well into effective matter dominated phase, any gravitational waves sourced via the 
turbulent mixing will be less strongly diluted. Consequently, the spectrum computed 
here represents a lower bound on the total signal. Further, the background will be 
strongly inhomogenous on small scales due to the bubble collisions and other preheating 
mechanisms. It is unclear exactly how efficient this preheating is and so how much it 
changes the effective equation of state, hence the last factor in Eq. (31) which may be 
significant. Thus, the suppression factor could again be reduced, leading to a stronger 
gravitational wave signal in the present epoch. These questions will most hkely be settled 
via direct numerical simulations, and we will tackle these projects in future work. 
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Appendix A. Thermal tunnelling of the flaton 



The bubble nucleation rate is [45, 46, 47] 



r ~ exp 



T 



where 



S,{T)^S[Ur),T]-S[0,T] 

poo 

S[(j){r),T] = 47r / r^dr 
Jo 



(A.l) 

(A.2) 
(A.3) 



is the energy of a nucleated bubble, and 0c is the critical bubble (bounce solution) which 
minimizes the energy, (f)^ satisfies 



dV 
90 







or0c _^ 2d0c 

dr"^ T dr 
with boundary conditions 

dr 

0c(oo) = 

Before percolation, the fraction of the universe consumed by bubbles is [49] 



and the time of percolation, tp, is given by 



(A.4) 

(A.5) 
(A.6) 

(A.7) 
(A.8) 
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During thermal inflation 

a{t) = e^<=* (A.9) 
where He is given by Eq. (17). Therefore Eq. (A. 7) becomes 

F(t) = % f dt'V{t') ^ [e^=(*-*') - ll' (A.IO) 



oo c 



Assuming + 



He HcT 



> 1 (A.ll) 



oo 
oo 

dss^r{t-s) (A.14) 







then Eq. (A.IO) is dominated by bubbles nucleated at times t' well within a Hubble time 
of the final time 

H,{t-t')<^l (A. 12) 

and so 

F(t) £ dt'mit-t'f (A.13) 

~ T 

where s = t — t'. Taylor expanding InF about the time of dominant bubble nucleation 
t' — tn gives 

T{t - s) = r(t,) exp [- /3{t^) (s-t + t^) + ...] (A.15) 
where the higher order terms are small if * 

|/?| (f - tn)' ~ 1^ « 1 (A. 16) 

Then Eq. (A.14) becomes 

F{t) =— dss^ r(tn) exp [- P{Q {s-t + 1^)] (A.17) 
<J Jo 

^ ^I(^e^(*n)(t-*n) r (A.18) 

3 P{tn} Jo 



where u = P{tn) s. Therefore 

87rr(t„) 



with the dominant contribution coming from u — Therefore the dominant bubbles 
at percolation are nucleated at 

(„ = tp-^ (A.20) 

+ Sec Appendix D for numerical estimates. 

* Numerical estimates give |/?|//3^ ~ 10~^ to 10~^. 
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and their size at percolation is 
R ' 



FYom Eqs. (A.8), (A. 19) and (A.20), i„ is given by 



r(tn) 



87re3 



(A.21) 



(A.22) 



Taking td = reexpressing the time in terms of the temperature and using Eq. (A.l), 
we define the /3 parameter used in the body of the paper as 



(A.23) 



/3 _ 


r 




^ d 






Ho 


HT 








Tci 



with Tci determined by solving 

SsiTd] 



exp 



cl 



87re3 V H, 



(A.24) 



Appendix B. Flaton finite temperature effective potential 



Before the phase transition, the flaton (p is strongly coupled to the thermal bath and its 
finite temperature effective potential is given by 

y(0,T) = V^(°)(0) + y«(0,T) 



(B.l) 
(B.2) 

(B.3) 



where the sum is over all the particles in the thermal bath, Qp is the number of degrees 
of freedom of the pth particle, Jp — J± for bosons and fermions respectively, and [50] 



(B.4) 



The thermal effective mass-squareds of the non-MSSM quark superfields (which 
dominate the fiaton's couplings to the thermal bath - see Appendix C) are [51, 52] 

where we assume the boson masses ~ rUg and the fermion masses are negligible. 
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Appendix C. Flaton couplings to the thermal bath 

The flaton must have unsuppressed interactions with the thermal bath in order to be 
held at the origin during thermal inflation. It is a gauge singlet so the interactions 
should be Yukawa couplings in the superpotential 

J]Ai$^i^, (C.l) 

i 

where $ is the flaton superfield. After thermal inflation, once the flaton has reached its 
vacuum value (f) — (f)o, the superfields and acquire masses 

M^^ = ^ ~ 10^° to 10^2 GeV (C.2) 

and so are not part of the MSSM. 

The and should form complete representations of SU(5) in order to preserve 
gauge coupling unification, and the size of the representations should be limited to at 
most two lO's and two lO's in order to preserve perturbative gauge coupling unification 
[53]. The and affect the renormalisation of the gauge couphngs and their 
masses are proportional to (p, therefore the gauge couplings will acquire a logarithmic 
dependence on 0. In particular, the strong gauge coupling during thermal inflation, 
g3{0), is related to the strong gauge coupling in the vacuum, g3{(f)o) — 1-2, by 

where Sriq is the number of colored superflelds in "if. For example, = 1 for = 5 
and = 3 for ^ = 10. 

At least some of the Yukawa couplings Aj should be strong enough to drive the 
fiaton's mass-squared negative at low energies. The renormalisation of the Yukawa 
couplings is given by 

dXj Ao 



din A 167r2 



SXl + J2x'j-^J2Ca(Ri)9l (C.4) 

where Ca{Ri) is the quadratic Casimir invariant for the superfield i in representation Ri, 
for example Ca = 4/3 for a superfield in a 3 of SU(3). Therefore the quark superfields 
will have the largest Yukawa couplings. However, quarks in different multiplets could 
have different Yukawa couplings, so the number of quarks with strong Yukawa couplings 
n\ could be less than the total number of quarks rig. Strong quark Yukawa couplings 
will be pulled close to their renormalization group infrared fixed point 

{Sn, + 2)Xl^ = -gl{0) (C.5) 

where Sux is the number of colored superfields with strong Yukawa couplings and App 
is the infrared fixed point of their Yukawa coupling. Therefore, from Eq. (C.3), 

App = '^^^'^"^ (C.6) 



^(9n, + 6) [l + ^^|(0o)ln 
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Thus the strongest Yukawa couphngs can be estimated and these will be the couplings 
that dominate the coupling of the flaton to the thermal bath. Values of App for the 
plausible choices of representations are shown in Table CI. 





5 


2x5 


10 


5 + 10 


2 X 10 




1 


2 


3 


4 


6 


93(0) 


1.05 


0.94 


0.85 


0.79 


0.69 




5 


5 


2x5 


10 


5 + 10 


10 


2 X 10 


nx 


1 


1 


2 


3 


4 


3 


6 


App 


1.09 


0.97 


0.77 


0.59 


0.49 


0.48 


0.36 



Table CI. The infrared fixed point Yukawa couplings, App, for (/)o = 10^^ GeV and 
nis = 10^ GeV. ^'efj are the representations with strong Yukawa couphngs. 

In the case of the baryogenesis scenario of Ref. [25], the Affleck-Dine field LHu has 
an expectation value /q ~ 10^ GeV at the end of thermal inflation, temporarily giving a 
large mass to some of the MSSM supcrflelds. Eq. (C.3) then becomes 

1 _ 1 1 
gl{0,lo) ~ 9l{cl>o,0)^ie^ 

with A„ ~ 10-^ Ac ~ 10-2 and A* ~ 1, and Table CI is replaced by Table C2. 
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5 + 10 


2 X 10 
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4 
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1.31 


1.11 


0.98 


0.88 


0.76 
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2x5 


10 


5 + 10 


10 


2 X 10 
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1 


2 


3 


4 


3 


6 


App 


1.35 


1.14 


0.90 


0.68 


0.54 


0.53 


0.39 



Table C2. The infrared fixed point Yukawa couplings, App, for (po = 10^^ GeV, 
nis = 10'^ GeV and Iq = 10^ GeV in the case of the baryogenesis scenario of Ref. [25]. 

Appendix D. Estimation of P/Hc 

We obtained the estimates for f3/Hc given in Table Dl by numerically integrating 
Eq. (A. 4) and solving Eq. (A. 23) for Td, using the flnitc temperature effective potential 
given in Appendix B with the gauge couplings g = g3{0) and Yukawa couplings A = App 
given in Table CI. 

In the case of the baryogenesis scenario of Ref. [25], using Table C2, Table Dl is 
replaced by Table D2. 
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0.5 


900 


1100 


1600 


3100 


5600 


6100 


16000 


1 


1000 


1300 


1900 


3800 


7800 


8300 


23000 


2 


1200 


1600 


2600 


5700 


12000 


14000 


39000 



Table Dl. Estimates of /?/i?c for Hc/Td = 10 is the ^ boson mass and we 

set the fermion mass equal to zero. 
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600 


800 


1100 


2000 


3700 


4500 


11000 
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700 


900 


1300 


2500 


4800 


5900 


17000 


2 


800 


1100 


1600 


3400 


7100 


9100 


26000 



Table D2. Estimates of (i/Hc in the case of the baryogenesis scenario of Ref. [25]. 
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